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Abstract — We develop a behavioral framework for 
switched linear differential systems (SLDSs). The definition 
of a SLDS is introduced, and several conditions for stability 
are derived using the notion of positive-realness. 

Index Terms — switched systems; behaviors; quadratic 
differential forms; positive-realness. 



I. Introduction 

In the classical approach (see for example J3], 0), 
switched systems consist of a bank of systems rep- 
resented in state-space or descriptor form (see 11181 ). 
together with a supervisory system that produces a 
switching signal determining which of the systems in 
the bank is currently active. When switching occurs, the 
continuity of the state trajectories may or not be required; 
in the latter case, "reset maps" are specified in order 
to produce new initial conditions for the post-switching 
evolution. 

In many physical situations, modeling switched sys- 
tems by means of state representations sharing the same 
state space is justified from first principles: for example, 
when dealing with switched electrical circuits, it is 
usually necessary to consider the state of the overall 
circuit in order to model the transitions between the 
different dynamical regimes. However, there are other 
situations in which modeling a switched system using 
representations that share a common, global state space 
is not justified from physical considerations; for example, 
when considering a control system consisting of a plant 
and a bank of controllers which may have different or- 
ders. In this case the dynamical regimes of the controlled 
system may have different degrees of freedom; it may 
be possible to conceive a global state space common 
to the different dynamics, but there is no compelling 
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reason to do so. A more fundamental issue with the 
classical approach is that modeling a physical system 
from first principles hardly ever results in a state-space 
description- this usually needs to be constructed from 
the set of higher-order differential equations describing 
the model. 

In this paper we propose a novel notion of switched 
system based on a behavioral point of view. In our 
approach a switched system consists of a set of dy- 
namical systems described by linear, constant-coefficient 
higher-order differential equations, together with some 
"gluing conditions" that relate the system variables and 
their derivatives immediately before and after switching. 
Thus, we do not assume that the different dynamics 
are represented in state form; moreover, the behaviors 
describing the different dynamical regimes are not even 
assumed to share the same state space; the minimal state 
space dimension for the systems may be different. 

A critical issue in switched systems is their stability, 
since switching among stable systems may give rise to 
instability. So many diverse solutions to this problem 
in the classical framework have been proposed that it 
is impossible to acknowledge all contributions in this 
area; we refer instead to ||6] for an exhaustive survey 
of the available results and for some novel ones. In this 
paper, we shall concentrate on the Lyapunov approach 
to stability: we show that the existence of a suitably 
defined higher-order Lyapunov function also ensures the 
stability of the switched behavioral system, and we study 
sufficient conditions for the existence of such Lyapunov 
function. 

The paper is organized as follows. The employed no- 
tation and a summary of the relevant background notions 
is given in section ITT] In section [III] we introduce the def- 
inition of switched linear differential structure (SLDS), 
and we discuss the concept of gluing conditions. In 
section [TV] we present a sufficient condition of Lyapunov 
type, based on the concept of strict positive-realness, for 
the stability of a SLDS consisting of two subsystems. 
In section [V] we introduce the notion of positive-real 
completion, and state a necessary and sufficient condition 
for the existence of a Lyapunov function for a SLDS with 
two subsystems. Some important remarks are provided 



2 



in section |VTJ The proofs of the most technical results 
are given in the Appendix. 

II. Background material 

A. Notation 

The space of n dimensional real vectors is denoted by 
R n , and that of mxn real matrices by R mxn . R mx ' denotes 
the space of real matrices with m rows and an unspecified 
finite number of columns, and R ,x * the space of real 
matrices with a finite but unspecified number of rows 
and columns. R mXO ° denotes the set of real matrices with 
m rows and an infinite number of columns. mX m denotes 
the zero m x m matrix, and I m denotes the m x m identity 
matrix. Given matrices A and B with the same number 
of columns, col(A, B) denotes the matrix obtained by 
stacking A over B. 

The ring of polynomials with real coefficients in 
the indeterminate £ is denoted by R[£]; the ring of 
two-variable polynomials with real coefficients in the 
indeterminates £ and i] is denoted by R[£, rf\. R rXw [£] 
denotes the set of all r x w matrices with entries in 
the ring R[£], and R nXm [C,r/] that of n x m polynomial 
matrices in £ and 77. The set of rational mxn matrices 
is denoted by R mXn (£). 

The set of infinitely differentiable functions from R 
to W is denoted by £°°(R,R W ). If / : R — > R*, we 
denote f(t~) := l\m T ^ t f(t) and f(t+) := lim rV /(t) 
provided that these limits exist. The set of linear dif- 
ferential systems with R w -valued signals is denoted by 
£ w . 

B. Linear differential behaviors 

A linear differential behavior is a linear subspace 
53 of £°°(R,R W ) consisting of all solutions w of a 
given system of linear constant-coefficient differential 
equations, represented by 

*(!)— 0, a) 

where R £ R ,Xw [£]. Equation © is called a kernel 
representation of the behavior 53 := {w G £°°(R,R W ) | 
w satisfies (Q~|)}, and w is called the external variable of 
53. The class of all such behaviors is denoted by £ w . 
Kernel representations are not unique: for example, if 
U € R* x, [£] is unimodular (see El for a definition) and 
53 = ker R (|), then also 53 = ker (UR) (f t ). 

When 53 is controllable (see Ch. 5 of ifTTIO . it can be 
represented in image form 

w = M (—] £ , (2) 



where M € R wX 1 [£] and £ is an auxiliary 
variable called the latent variable; i.e., 53 := 
{w I 3 £ such that © holds}. 

Given 53 € £ w , it may be possible to choose some 
components of the external variable w freely (in the 
sense of Def. 3.3.1 of ifTTTO : such free variables are then 
called input variables. It can be shown that the maximal 
number of input variables is an invariant, denoted by 
m(53). Once m(53) free variables have been chose, the 
remaining components of w are output variables; evi- 
dently, the number p(53) := w— m(53) of output variables 
is also an invariant. If 53 has no input variables it is called 
autonomous; in this case m(53) = 0, or equivalently 
p(53) = w, and it can be shown that 53 admits a 
kernel representation with R square and nonsingular. 
If 53 is controllable, by permuting the components of 
w an image representation © can be achieved with 
M = co\(U,Y) and U € R lxl [£] nonsingular; then 
u = U (4A £ is an input variable, and y = (-£) I is an 
output variable. Note that m(53) = 1, the dimension of £. 
The transfer function from u to y is the matrix of rational 
functions YU^ 1 . It can be shown that an input-output 
partition w = col(n, y) can be found so that YU^ 1 is 
proper. 

C. State maps 

A latent variable I is a state variable for 53 if and 
only if there exist E, F £ R ,x ', G G R ,Xw such that 
53 = {w I 3 £ s.t. Ef t + F£ + Gw = 0}, i.e. if 53 has 
a representation of first order in £ and zeroth order in 
w. The minimal number of state variables necessary 
to represent 53 in this way is an invariant called the 
McMillan degree of 53, denoted by n(53). 

In |fT3l it has been shown that a state variable for 
53 can be computed as the image of a polynomial 
differential operator called a state map. We now review 
the construction of state maps for autonomous behaviors. 

1) State maps for autonomous systems: In this case, 
a state map acts on the external variable w. Let 53 
represented as in CO with R € R WXw [£] nonsingular, and 
consider the set defined by 

X(R) := {/ € R lxw [£] I fR" 1 is strictly proper} . 

(3) 

3L(R) is a finite-dimensional subspace of the vector space 
R lxw [£] over R. To construct a state map X € R* xw [£] 
for 53, choose a set of generators x% € R lxw [£], i = 
1,...,N of 3L(R), and define X := col(xi)i=i ) ... ) jv', to 
obtain a minimal state map, choose {a;i}i=i,...,jv so that 
they form a basis of X(R). 

It can be shown (see Cor. 6.7 of |[T3l ) that the 
McMillan degree n(53), i.e. the minimal dimension of 
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a state variable for 53, equals deg(det(i?)). Note that if 
X £ IR ,Xw [£] is a state map for 53 = ker R (^), then for 
every F € R' Xw [f], also X'(£) := X(£) + F(£)R(£) is 
a state map, since X' (^) w = X (^) w for all td £ 53. 

Finally, given an autonomous behavior ker J? (-4) 
associated with i? nonsingular, it can be shown that there 
exists a state map X and A £ M***, F € M ,Xw such that 
£X(£) = + 

2) Functions of the state: Let 55 6 £ w , and let X £ 
IR ,Xw [£] be a state map for 53 acting on w. A polynomial 
differential operator d(4A is a (linear) function of the 
state of 53 if there exists a constant vector / £ R lxw 
such that d(i)w = fX H) w. 

3) R-canonical representatives: We recall the notion 
of ^-canonical representative of a differential operator. 
Let R £ M wXw [£] be nonsingular, and let / £ M lxw [£]; 
/ can be uniquely written as fR~ x = s + n, where 
s € M lxw (£) is a vector of strictly proper rational 
functions, and n £ M lxw [£]. We call sR £ M lxw [£] 
the (uniquely defined) canonical representative of f 
modulo R, denoted / mod R. Note that the polynomial 
differential operators / (^) and /'(^), with f = f 
mod R, are equivalent along ker R f^§) in the sense 

/(£)« = /'(I)™ for a11 w € ker R (m)- 

The definition of i?-canonical representative extends in 
a natural way to polynomial matrices. 



D. Quadratic differential forms and dissipativity 

The parametrization of quadratic functionals of the 
system variables and their derivatives using two-variable 
polynomial matrices has been studied in EOl . to which 
we refer for a thorough exposition. 

1) Quadratic differential forms: Let $ £ IR wXw [£,?/]; 
then $(C,r?) = £ M * M CV. where *M e M wXw 
and the sum extends over a finite set of nonnegative 
indices. 77) induces the quadratic differential form 
(QDF) acting on C^-trajectories defined by Q$(w) := 
T,h,k(^dw) T ®h,k^-- Without loss of generality one can 
assume a QDF to be induced by a symmetric two- 
variable polynomial matrix $(£,77), i.e. one such that 
^(C) 7 /) = ^( r l^C) T > we denote the set of such matrices 
by W 8 *"[C,Tj\- 

$(£,77) can be identified with its coefficient matrix 



h ,fc=0,...,oo 



, in the sense that 



d>(c,r/) = [j w a» 



Although <l> is infinite, only a finite number of its 
entries are nonzero, since the highest power of ( and 



77 in $(£,77) is finite. Note that is symmet- 

ric if and only if <I> T = <3>. In that case, $ can 
be factored as <l := M T S$M with M an infinite 
matrix with an infinite number of rows and all but 
a finite number of elements equal to zero, and a 
signature matrix. The latter decomposition leads after 
premultiplication by [/„ • • •] and postmultiplica- 

tion by col [J„ 77l w • • •] to the factorization 77) = 
M (C) T S<j,M(?7) which, if we take M surjective, is called 
a canonical factorization of $(£,77). 

Given Qif, its derivative is the QDF Q$ defined by 
Q$(w) := f t (Q^{w)) for all w £ er°°(M,M w ); this holds 
if and only if $((,r)) = (C + »?)*(C>»7) ( see El, P- 
1710). 

A QDF Q$ is nonnegative along 53 £ £ w , denoted 

58 

by Q$ > if Q$>(w) > for all tx> £ 53; and positive 

23 23 

a/o«g 53, denoted by > 0, if Q$ > and [Q${w) = 
;w G 23] ==> [u, = 0]. If 53 = C°°(M,R H ), then 
we call (5$ simply nonnegative, respectively positive. 
Characterizations of nonnegativity and positivity in terms 
of algebraic properties of the matrix $(£,77) are given 
in (201, pp. 1712-1713. 

An analogous concept of i?-canonical representative 
is introduced for two-variable polynomial matrices; two 
QDFs Q$ l , Q® 2 are equivalent along ker R{j%) if 
Q$ 1 (w) = Q<$> 2 (w) for all w £ ker R{4j)- Among 
all QDFs that are equivalent to a given Q$, there is 
exactly one which is /^-canonical. This canonical repre- 
sentative can be computed as follows. Write $(£,77) = 
M(C,) T N(rj) (such a polynomial factorization can be 
computed factorizing the coefficient matrix of $(£,77) 
as $ = M N). Now compute the /^-canonical repre- 
sentatives M' of M and N' of N. Then the i?-canonical 
representative of <!>(£, 77) is M' '(C) 1 ~N' '(77). 

2) Dissipativity and Positive Realness: A control- 
lable behavior 53 £ £ w is dissipative with respect to 
the supply rate Q§ if there exists a QDF Qq, called 
a storage function, such that Q®(w) — ^Q^(w) > 
for all 77; £ 53. This inequality holds if and only if 

23 

there exists a dissipation function, i.e. a QDF Q/\ > 
such that for all w £ 53 of compact support it holds that 
f^Qz(w)(t)dt = /+~ Q A (w)(t)dt (see Prop. 5.4 
of |[20l ). Moreover, there is a one-one correspondence 
between storage- and dissipation functions, defined by 
f t Qq,(w) + Qa{w) = Q<s>{w) for all 77; £ 53. If 
53 = £°°(IR,IR W ), this equality holds if and only if 

(C + 77)*(C,7?) + A(C,77) = <I>(C,r?). 

Recall the notion of strict positive-real rational func- 
tion (see Q): a strictly proper rational matrix G := 
R2R1 1 , with Ri,R 2 £ M wXw [£], is strictly positive real 
if G(£) is real when £ is real, G is analytic in the right- 
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half plane, and G(-joj) T + G(ju) > V u G R. 

III. Switched Linear Differential Systems 

A switching linear differential structure is a finite 
family of linear differential systems, together with a set 
of switching signals that specify which system is active 
at each time instant, and a set of "gluing conditions" 
regulating the concatenation of the trajectories at the 
switching instants. 

Definition 1. A switching linear differential struc- 
ture (SLDS) £ is a quadruple £ = {V,J-,S,Q} 
where V = {1,...,N} C N, is the set of indices; 
T = (<8 1 ,...,«8 JV ), with <Bj 6 £ w for j G V, 
is the bank of behaviors; S = {s : R — > V 
s is piecewise constant and right-continuous}, is the set 
of admissible switching signals; and 

g (R' xw [e]) 2 } , 

with (k,£) £ V x V, k ^ £ is the set of gluing 
conditions. The set of switching instants associated with 
s G S is defined by T s := {t G R | lim r ^ s(t) ^ 
s(t)} = {h,t 2 , ■ ■ ■}, where U < t i+1 . 

A SLDS induces a switched behavior, defined as 
follows. 

Definition 2. Let S = {V, F, S, Q} be a SLDS, and let 
s G S. The s-switched linear differential behavior 93 s 

is the set of trajectories w : R — > R w that satisfy the 
following two conditions: 

1) for all ti,ti + \ G Y s , there exists 93 ^ G T, k G V 

such that H[ti,t i+ i) e ®fc I [**,**+!)»" 

2) to satisfies the gluing conditions Q at the switching 
instants for each ti G T s , i.e. 



The switched linear differential behavior (SLDB) 93 s of 

S itf de^nerf by Q3 E := U se5 

In the rest of this paper we make the standard as- 
sumption (see e.g. sect. 1.3.3 of |fT71 ) that the number 
of switching instants in a given finite interval of R is 
finite. 

The polynomial matrices GjJ" f (£), G^ £ (^) in Def.s 
Q] and [2] specify concatenation conditions on w at the 
switching instants. These conditions are an essential 
component of the dynamical laws of a switched linear 
system to be analyzed and are defined according to the 
particular nature of the switched system under study. 



Example 1. Assume that depending on the value of 
the switching signal, a plant Sp with two external 
variables, described by the differential equation iwi — 
w\ — W2 = 0, is connected with one of two possi- 
ble controllers T,c 1 and Sc 2 , described respectively by 
and — 2w\— u>2 = 0. Depending 



on which controller is active, the resulting closed-loop 

"I- 1 ~] 
1 d 



behaviors are 25 1 
1 



ker 



o d 



and 2? 2 := 



kcr 



dt 



. Note that 23 1 and 2^2 have different 



dt x dt 

1 

2 -1 

McMillan degree (2 and 1, respectively). We define the 
gluing conditions for the SLDS associated with 93 1 and 
"0 1 " 

-2 

and by G^ 2 (0 := [6 l], G+^ 2 (0 := [0 



93 2 by 



1 

1 
II. The 



rationale underlying our choice of gluing conditions is 
that any trajectory of 93 1 is uniquely specified by the 
instantaneous values of col(wi, W2), while a trajectory 
of 93 2 is uniquely specified by the instantaneous value 
of W2- Thus, when switching from the dynamics of 93 1 to 
those of 932, we require that the values of W2 before and 
after the switching instant coincide. In a switch from 232 
to 93 1, since the second differential equation describing 
232 yields W2 = —2w± before the switch, we impose that 
w 2{t~k) = W2{t~ k ) = — 2wx(tT). Note that this makes 
the switched trajectory as smooth as possible, taking 
into account the restrictions imposed by each individual 
behavior 23 1 and 932- □ 

While in principle Def . [Qdoes not restrict the choice of 
gluing conditions, in practice they should be well-posed, 
according to the following definition. 

Definition 3. Let S be a SLDS as in Def UI Let G 

R* Xw [£] be a minimal state map for 23 j, i = 1,...,N. 
The gluing conditions Q are well-posed if for all k, I = 
1, . . . , N with k ^ I there exists a map Lk^e : R* — > R* 
such that for all ti G T s 

[s(U-i) = k,s(ti)=£] 
d 



and 



d 



Thus, if G is well-posed then after each switching 
instant the state of the active behavior is specified as a 
function of the state before the switch. This property ex- 
presses the reasonable requirements that after the switch: 
a) it should be possible to uniquely define (together with 
the free components of w) the future evolution of a past 
trajectory; and b) that there be a causal link between the 
future and the past evolution. 
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Example 2. The gluing conditions of Example Q] are 
well-posed. A minimal state map for 23 1 is Xi(^) := I2, 
and a minimal state map for 232 is -X^CO = [l 0]. 
Def. [3] is satisfied with L\^2 '■= [l 0] and £2^1 := 
[1 2] T . □ 

In the rest of this paper we focus on the stability 
problem for SLDS, and we consider SLDS consisting of 
two autonomous behaviors: thus in the following F = 
(231,532), with 23^- = ker Rj (|), where Rj G R wXw [£], 
j = 1, 2, are nonsingular polynomial matrices. Moreover, 
we assume that Ri is column reduced (see (H), and 
that R2R1 1 is strictly proper. The first assumption is not 
restrictive, since a column-reduced kernel representation 
of 23 1 (for example, the row-Hermite form described in 
can be obtained by premultiplication of an arbitrary 
kernel representation by a suitable unimodular matrix. 
The second assumption has important consequences; to 
state them we need first to introduce the notion of stan- 
dard state map. Let 23 = ker R (£) with R G R wXw [£] 
nonsingular, and assume that R is column-reduced. De- 
note the column degrees of R by Sj, j = 1, . . . , w, and 
consider the polynomial vectors 

-k 



1, 



, w , k = 0, . . . , Sj — 1 



(4) 



where ej is the j-th vector of the canonical basis of 
R lxw . It follows from the material in section HT-Cll that 
the vectors dD form a basis for the state space of 23. 
Any matrix whose rows are the vectors (@) (or any 
permutation thereof) is called a standard state map for 
23. The following is a crucial result used repeatedly in 
the rest of the paper. 



Lemma 1. Let 23j = ker Ri (^), i = 1,2. Assume that 
i?i,i?2 G R WXw [£] are nonsingular, that R± is column 
reduced, and that R2R1 1 is strictly proper. Denote 
Hi := deg(det(i?j)); then n 2 < m, and there exist 
X[ G R( ni - n2 ) xw [£], X 2 G R naXw [£] such that X 2 (£) is 
a standard state map for 232, cmd the polynomial matrix 
Xi(-j-) defined by 



*i(0 == 



(5) 



is a standard state map for 23 1. Moreover, there exists 

n e M(«i-«2)xn 2 such that x ^ mod R2 _ ux 2 {0- 

Proof: The inequality 112 < n\ follows straight- 
forwardly from the following result and the material in 
section III-Cll 



Proof: See the Appendix. ■ 

To prove the claim on X\ defined by ©, bring 
i?2 via left unimodular multiplication to a column 
reduced form R' 2 (for instance, into row-Hermite 
form). Construct a standard basis for the space 
j/ G M lxw [^] j fR'2 1 is strictly proper |, and arrange 

the vectors of this basis in a matrix -X"2(£); then X2 (i) 
is a standard state map for 232 = ker R' 2 (4A = 
ker R2 (j?). Since R\ is column reduced and R2R1 1 
is strictly proper, applying Lemma |2] again we conclude 
that the rows of X2 (^) are a subset of the rows of any 
standard state map X\ (f t ) for 23i. Now rearrange the 
rows of Xi(£). 

To prove the existence of IT, observe that each row 
of X[ mod i?2 belongs to the state space X2 of 232 = 
ker i?2 (4f)> an d consequently can be written as a linear 
combination of the rows of Xi- This concludes the proof. 



Example 3. When w = 1, R\ is trivially column- 
reduced, and R2R1 1 is strictly proper if and only if 

n\ = deg(i?i) > deg(i?2) = "-2- A standard state map 
for 23 1 = ker R\ (4) consists of £ k , k = 0, . . . , n\ - 1, 
whose first ri2 elements form a standard basis for the 
state space of 232 = ker R 2 (4) . The rows of the matrix 
IT defined in Lemma Q] consist of the coefficients of the 
polynomials £ fc mod i?2(£)> ^ = n 2> • • • > n\ — 1. □ 

In the rest of this paper we consider standard SLDS, 
defined as follows. 



Definition 4. Let £ = {P,F,S,G} be a SLDS with 
T = (23i, 23 2 ), where 23j = ker Rj (f t ), j = 1, 2, with 
Rj £ M wXw [^] nonsingular polynomial matrices, R\ col- 



umn reduced, and R2R1 1 strictly proper. Denote 



n : 



deg(det(-Rj)), j = 1,2, and let X[ G R(»i-n») x »[f], 
X 2 G R n2Xw [£] and n G £ e de/ined as m 

Lemma \T\ £ is a standard SLDS if ffte gluing conditions 



are 



and (G 



1(6) 

-*(0) : = 



^2(6 
nx 2 (o 

(* 2 (0, *>(£)). 



^2(0 



It is a matter of straightforward verification to check 
that the gluing conditions of a standard SLDS are well- 
posed. 



Lemma 2. Let i?i,i? 2 G 



i?2^i 1 is strictly proper, and define Xi := {/ G 
R lxw [£] I fRi 1 is strictly proper}, i = 1,2. 77zen X 2 C 



'[^], such that Example 4. Consider the case of Example |3l and assume 
that Ri and R2 are monic, and that n\ = n 2 + 1. Denote 
^2(0 =: EjLo 1 R2,j&- The gluing conditions of 
the corresponding SLDS are (G^iCO) ^2^1 (0) = 



6 



V 



-R 



2,0 



£fti-l 
" 1 

n-L-2 



and 



ni-2 



When switching from 252 to 25i, to obtain "initial 
conditions" uniquely specifying w G 25 1, we need 
to define the {n\ — l)-th derivative of w after the 
switching instant. The gluing conditions impose that 
this coincides with the value of the (rti — l)-th 
derivative of w before the switching, given by the 
equality ^k^w = -R 2 ,ow - ... - R^n-2§0iW 
arising from the fact that w G 252 before the switch. 
When switching from 25 1 to 252, we project the vector 
of derivatives characteristic of a trajectory w G 25i 
down onto the shorter vector of derivatives characteristic 
of a trajectory w G 252- Note that our choice of gluing 
conditions guarantees that the corresponding trajectory 
is as smooth as possible. □ 

IV. POSITIVE-REALNESS AND LYAPUNOV STABILITY 
OF SLDS 

We call a SLDB 25 s (and by extension, the SLDS 
E) asymptotically stable if lmi£_ s . 0O w(t) = for all 
w G 25 s ; note that if 25 s is stable, all the 25, in F 
are autonomous, since they cannot have free variables. 
Our notion of stability, under arbitrary switching signals , 
is analogous to the notion of uniform stability of the 
classical switching systems literature, see section 2.1.1 
of 0. 

When dealing with linear differential behaviors, 
asymptotic stability can be proved by producing a 
higher-order Lyapunov function (see 111011 . 11201 ). i.e. a 

23 

quadratic differential form Q§ such that Q$, > and 

, 03 

;|Q$ < 0, see Th. 4.3 of GOl . In the classical framework 
for switched systems it is often difficult to find a common 
Lyapunov function, and multiple Lyapunov functions are 
used instead (see section III.B of (6l). This situation 
arises naturally in our framework, where the McMillan 
degrees of the behaviors of the SLDS are in general not 
equal; hence the following definition. 

Definition 5. Let S be a SLDS as in Def [7] A QDF 
is a Lyapunov function for £ if is a Lyapunov 
function for 25j = ker Ri(^), where ?y) = 

Vl/(Cj r/) mod Ri, i = 1, . . . , N; and moreover, for every 
w G 25 s and every switching instant it holds that 

Q*.^,) («>)(**) >Q* s(tfc) H(t+). 



The crucial requirement in Def. |5] is that at each 
switching instant the value of Qq, does not increase; it 
follows from standard arguments (see Th. 4.1 of ED) 
that if a SLDS £ has a Lyapunov function, then £ is 
asymptotically stable. 

This section is devoted to the proof of a sufficient 
condition for the existence of a Lyapunov function for 
a standard SLDS. In the following result, we associate 
to the matrices R\ , R2 in Def. [4] the transfer function 
R2R1 1 of the system described in image form by 
col(i?i, R2). The next result states that if R2R1 1 is 
strictly positive -real, then there exists a storage function 
Qq,, induced by a symmetric two-variable polynomial 
matrix ^ G M^ Xw [C,r?], associated to R2R1 1 that is a 
Lyapunov function for the SLDS. 



Theorem 1. Let 

{25i = ker R x (|) , 25 2 = ker R 2 (|) }. Defi 



£ be a standard SLDS with T 
= Ker rt2 1 ejine 



Tl 



W / w 






I^j W 







= R^C) 1 R 2 {rj) + R 2 (C) J Ri(v) ■ 

Assume R2R1 1 is strictly proper and strictly positive- 
real. There exists A € M wXw [^] nonsingular such that 
det(.A) is an anti-Hurwitz polynomial, and AR^ 1 is 
strictly proper. Moreover, 

^(C^-AiCVAir,) 



*(C,r?) :-- 



eClCr)], (6) 



induces a Lyapunov function for S. 

Proof: R2R1 1 being strictly positive real implies 
that $>(—ju,joj) > for all oj G M. By standard 
results on polynomial spectral factorization (see |f2]), we 
conclude that there exists A G IR wXw [£] such that det(A) 
is anti-Hurwitz and -£) - A(-£) T A(£) = 0. This 
implies that rj) — A(() T A(r]) is divisible by (C + v) 
and therefore there exists ^((,7]) as in Q. We prove 
that is a Lyapunov function for E by showing that 
it satisfies the conditions 

23! 

CI. > and f t Q^ < 0; 

C2. Qy > and ^Q* < 0; 

C3. Qq, does not increase at switching instants. 

To prove CI, first use the fact that strict positive -realness 

implies that the largest storage function (corresponding 

to the anti-Hurwitz factorization of $(—£,£), see Th. 

5.7 of EOlO is nonnegative definite, i.e. Qq, > (see 

. . 23i 

11201 . Theorem 6.4); this implies a fortiori that Qq, > 0. 
To prove the second claim of CI, note that from © it 



follows that j^Qm 



* A (|) .This functional 



is strictly negative along 25 1, since A is anti-Hurwitz and 
Ri is Hurwitz. 
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To prove that condition C2 is satisfied, we first com- 
pute the canonical representative of 77) modulo i?2, 
*'(C,r?) := *(C^) mod #2- Since = Q*, and 

Q* > 0, it follows that Q^> > 0. To prove the second 
part of C2, observe that the canonical representative of 
the derivative of Qq,> is induced by (( + ^^'(Cr/) := 
{( + V)y((,v) mod R2 = -A'(Q T A'(r]), where A' 
is the canonical representative of A modulo R 2 . Con- 

sequently, to prove that j^Q^' < it suffices to show 
that col(A'(A), i? 2 (A)) is full column rank for all A G C, 
which guarantees that ^ (Q<ij>(w)) is non zero for the 
trajectories of 232- Assume by contradiction that there 
exists A* G C such that col^^A*), R 2 (X*)) is not full 
column rank. Then i?2(A*) is not full rank and, since 
i? 2 is Hurwitz, we conclude that A* G C~. Also A'(X*) 
is not full rank, and since A(£) = A'(£) + P(£)i? 2 (£) 
for some P G M wXw [£], this implies that A(\*) is not 
full rank. Since det(A) is anti-Hurwitz, this yields a 
contradiction. Condition C2 is satisfied. 

We now prove that condition C3 is satisfied; the 
following fundamental fact will be instrumental to this 
purpose, and will be repeatedly used in this paper. 



This inequality can be rewritten as 



[x[{i) w {r k )\ 

Since X 1 (f t )w(t^) 



l 7l2 









L n 1 

n 



(V) 



> 



ni due to the minimality of X\{^ 



*11 L 



is equivalent with ^ 
Now partition \I> = 



is arbitrary in 
inequality (0 

^n 2 > 



*i" 2 



pri2 xrt2 



'r»2 



*12 
^22 

*12 G R njX ( ni - Bj ) and f 22 G I 
standard linear algebra manipulations, we obtain that (JVj 
is equivalent with 



n 

, with G 

(n 1 -n 2 )x(n 1 -n 2 ) After 



-(^12 + n' ^22)^22 (*i2 + ^2 2 n) 
^22 



> 



A standard linear algebra argument shows that the in- 
equality holds if and only if fyj 2 + ^! 22 Il = 0. In order 
to prove that this is the case, consider first the following 
lemma. 



Lemma 4. Define k 



then k G 



(tii-n 2 )xw 



Lemma 3. Let N,D G M wXw [£] with D nonsingular. Let 
ND^ 1 be the transfer function of the system described 
in image form by co\{D, N). Assume that ND^ 1 is 
strictly proper and strictly positive real. Define the k T (^>J 2 X 2 (£ i ) + ^f 22 X[ (£)) 
supply rate := D(() T N(r]) + N(() T D(i]). Let 

A(C, 7]) be a corresponding dissipation function, and 
write A(C,r/) = Q{() T Q{r]), far some Q G R ,Xw [£]. 
Then QD~ l is strictly proper. Moreover, let Qq> be a 
storage function, and let ^(C, 77) = M(^) T Sij,Af (77) 
be a canonical factorization of ^(^,77), for some M G 
IR* Xw [£] and some signature matrix Then MD^ 1 is 
strictly proper. 



Moreover, R 2 (£) = 



Proof: See the Appendix. ■ 
It follows from Lemma |4] and the last statement of 
Lemma Q] that 

= i? 2 (6 mod R 2 



k T (*J 2 X 2 (£) + * 2 2*i(0) mod ft 



(8) 



k ^ 2 + * 22 n )x 2 (0 . 



Proof: See the Appendix. ■ 

It follows from the material gathered in sections III-C 1 1 
and III-C2I and from Lemma [3] that having chosen a 
standard state map for 23 1 according to Lemma Q] 
we can_ write ^((,rj) = Xi(Q T ^Xi(r]) for some 



* 1 G 



ini xrti 



. Finally, since AR 1 1 is strictly 
proper, there exists A = A T G K niXni such that 
A(C,r?) = A(C) T A( V ) = X 1 (C) T AX 1 ( V )- 

To prove that C3 is satisfied first consider the case 
of a switch from 23 1 to 03 2 at t k . Taking the standard 
gluing conditions into account, we want to prove that 



[X[{i)w{r k ) 



nx 2 (|Mt+). 



> 



Now recall that the rows of X 2 (£) are linearly indepen- 
dent over M, since they are part of a minimal state map. 
Consequently ([8]> implies k T (^J 2 + ^ 22 U) = 0. We now 
show that k T G K wX ( ni_n2 ) i s square and nonsingular; 
this will lead to the desired equality ^2 + ^22!! = 0. 

Lemma 5. Define k as in Lemma @ Then n\ — n 2 = w, 
and k is nonsingular. 

Proof: See the Appendix. ■ 
From Lemma [5] we conclude that the inequality © is 
always satisfied, and consequently does not increase 
when switching from 23 1 to 232- 

To conclude the proof of the Theorem we need to 
show that the value of Qq, does not increase when we 
switch from 232 to 23i. It is a matter of straightforward 
verification that for this type of switch, the standard 



s 



gluing conditions imply that the value of Q$(w) before 
and after the switch is the same. The claim of the 
Theorem is proved. ■ 

Remark 1. QDFs act on £°° -functions, while the tra- 
jectories of a SLDS are in general non-differentiable; 
however, this mismatch in the degree of differentiability 
is relevant neither to Th. [T] nor to the other results 
presented in this paper, since we only use the calculus of 
two-variable polynomial matrices associated with QDFs 
as an algebraic tool. For example, in the proof of Th. Q] 
when considering the value of the functional Qq, before 
and after a switch, only the properties of its coefficient 
matrix are being used. □ 

Remark 2. Theorem Q] reduces the problem of finding 
a Lyapunov function for a standard SLDS to performing 
a polynomial spectral factorization, for which standard 
procedures can be used, and computing ^(£,77) from 
formula ©. □ 

Remark 3. In the scalar case (w = 1) the proof of Th. 
Q] simplifies considerably; from standard results in the 
theory of scalar rational positive-real functions it follows 
that the difference in the (McMillan) degrees of Ri and 
R2 is one, and several of the intermediate steps necessary 
in the proof of the multivariable case are trivial. See [12 1 
for details. □ 

Remark 4. Theorem Q] holds true also if the assumption 
on the strict properness R2R^ 1 is relaxed to i^-R^ 1 
being bi-proper, i.e. proper and with a proper inverse. 
It follows from the material in section III-C1I that in 
this case the state space of ker Ri(4f) and that of 
ker i?2 (4f) are the same. Choose the gluing con- 
ditions to be (^(0.^2(0) : = (X(0,X(0l 
{G^iO^LiiO) ■= (X(£),X(e)). It is a matter 
of straightforward verification to check that the largest 
storage function for R^R^ yields the desired Lyapunov 
function. For the scalar case (w = 1), this result has been 
proved in OH. □ 

Remark 5. In the state-space framework for switched 
systems it is well-known that if the open-loop transfer 
function of a system is positive-real, then all closed- 
loop systems obtained from it by state feedback share a 
quadratic common Lyapunov function (see section 2.3.2 
of @ and lfl31l . lfl6l ). Th. Q] offers a new perspective 
on the relation between positive-realness and stability 
of switched systems: in our framework, the different 
dynamical regimes do not arise from closing the loop 
around some fixed plant, and positive-realness arises 
from the interplay of the dynamics of the different 
systems. □ 



Remark 6. Positive realness is only a sufficient con- 
dition for the stability of switched differential systems, 
and not surprisingly it is not necessary. Define R\{£) := 



2 





£ + 2 



£ + 3 
-5 





e + i 
9 
4 



It is readily 



w 2 \ 



is a 



verified that V{w\,W2) '■= [wi 

common Lyapunov function (in the classica 
ker Riijg), i = 1,2; however, R2R1 1 is not positive 
real. □ 



sense) for 



V. Positive Real Completions 

Given a standard SLDS, the strict positive-realness 
assumption on it^^f 1 is in general not satisfied. The 
idea we pursue in this section is to "complete" R2 
premultiplying it by M € M wXw [£] to obtain a strictly 
positive-real transfer function {MR2)R[ ; M will be 
called a positive real completion of R2. If such M 
exists, then Th.Q] guarantees that there exists a Lyapunov 
function for the SLDS E' with behaviors ker i?i (^) 
and ker (MR2) (4z)\ in this section we show that a 
Lyapunov function exists also for the SLDS E with 
behaviors ker Ri (^A and ker R2 (4f)- Moreover, we 
also show the converse: if a Lyapunov function for E 
exists, then there also exists a positive-real completion 
for R 2 . 

Definition 6. Let T,bea standard SLDS. M € M wXw [£] is 
a strictly positive-real completion of R2RY 1 if MR2RY 1 
is strictly positive-real. 

Remark 7. The existence of a positive -real completion 
for given i?i,i?2 £ M wXw [£] can be checked computa- 
tionally solving an LMI associated with the dissipation 
inequality of the M parameter-dependent system repre- 
sented in image form by col(Mi?2, Ri)- Such problem 
is straightforwardly solved in the scalar case, where an 
upper bound on the degree of the polynomial M € 
namely deg(M) < deg(i?i) — deg(i?2)> is known. In 
the multivariable case unfortunately an upper bound on 
the degree of M is not available, and the need arises 
to iteratively solve LMIs of increasing dimensions. Note 
that such iterative procedures are standard when dealing 
with higher-order dissipative systems, see e.g. p. 32 and 
Algorithm 4.1 p. 39 of OH. □ 

The main result of this section is the following. 

Theorem 2. Let E be a standard SLDS. Assume ri2 < 
n\ — w. The following statements are equivalent: 

1. There exists a Lyapunov function Qq, for E; 

2. There exists a strictly positive-real completion M £ 
W**[£]for R2R1 1 . 
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Proof: ((1) ==> (2)) Without loss of generality, we 
assume Qq, in canonical form modulo R\. To compute 
a positive-real completion, we need first to decompose 
Qq, in a convenient way. Use Lemma Q] to define a 
state map X\ (4) for 05 1 as in ([5]), and Lemma [3] to 



conclude that there exists ^ 
that ^(C,rj 



such 



*7 2 



*22 

and ^22 € 



= Xi(C) T ^^i(r?). Partition * as * =: 
where #n € M n2Xn2 , € R n 3x(m-n 2 ) 



ri2 > . Consider the gluing con- 
ditions of the SLDS; following the same argument used 
in showing condition C3 in the proof of Theorem Q] 
we conclude that Qq, does not increase at the switching 
instants if and only if tyj 2 + ^22^ = 0. Consequently, 
we can rewrite >f as 



$ = 






-11*22 


-v 22 n 


^22 




0" 




" n T 




o_ 


+ 





-n 2 )_ 



^22 [n 



'(n 1 -n 2 )\ j 



(9) 



where we define vp' := — H T ^22^- Pre- and post- 
multiply © by X 1 (C) T on the left and Xi(rj) on the 
right to obtain 

*(C,r/) = 

X 2 (C) T ^'X 2 (r ? ) 



+ *i(0" 



-I 



n 

(ni-n 2 ) 



IT satisfies X[ (£) mod i? 2 
consequently 



*22 [n -I {n ^ n2) ]X 1 {r } ). 

(10) 

ILY 2 (£) (see Lemma [1}; 



[n 



'(rc-i— n 2 )J 



*2 (f ) 



w 







for all w £ 03 2 . We conclude that Qii- = Q^s and 
following an argument analogous to that used to prove 
C2 in Theorem [Q that the functional Qq,> induced by 
the first addendum *'(C,7/) := X 2 {() T &X 2 (ri) of (H0]» 
is a Lyapunov function for 05 2 . 

Having decomposed *(£,t/), we aim at rewriting the 
expression in order to make evident what the positive- 
real completion matrix M is. To this purpose, consider 
the following result. 

Lemma 6. Define 

P(0 := - (X[(Q mod R 2 ))R2(ar 1 

Then 

f (C,77) = *'(C,»7) + #2(0^(0^22^(^2(7?) • 

(11) 



Proof: See the Appendix. ■ 
We use (fTTT) to obtain a convenient expression for 
the two-variable polynomial matrix corresponding to the 
derivative of Qq, along 05 1. 

Lemma 7. Let P be as in Lemma [6] There ex- 
ists B € M("i-" 2 ) xw smc/i that -A(C,??) := (C + 
r/)^(C, rj) mod i?i satisfies 

-A(C, r?) =(C + »7)*(C, ^) " Ri(0 T B T ^22P(r j )R 2 (v) 
-R2(C) T P(v) T ^22BR 1 (rj) . 

Proof: See the Appendix. ■ 
From the proof of Lemma [7] it follows that B := 

lim^ oo ^P(0-R2(0#i(0~ 1 - We Ciarm that the matrix 
M defined by 



M(0 := S T * 22 P(0 , 



(12) 



is the desired positive-real completion of R 2 R 1 1 . Indeed, 
the equality established in Lemma |7J together with 

23 1 

— A(C 77) < and i?i being Hurwitz, prove strict 
positive -realness of MR 2 R~[ l . Moreover, MR 2 R^[ 1 is 
strictly proper, since 

b t * 22 p(0R2(0Ri(0~ 1 

= B T * 22 (X{(£) - (X[(0modR 2 ) R^O^MO) MO' 1 
= B T ^ 22 (X[(0 - (X[(0 mod R^R^r 1 ) , 

and since X[(£) - (X[ (0 mod #2)^1 (f)" 1 is strictly 
proper. This concludes the proof of ((1) =^> (2)). 

((2) =^> (1)J. To prove the claim we show that 
there exists e € M, e > such that among the storage 
functions for (eMR^R^ 1 there is one satisfying the 
conditions CI, C2 and C3 stated in the proof of Theo- 
rem [T] Such storage function also induces a Lyapunov 
function for S. 

Lemma 8. Let (MR^R^ 1 be strictly proper and strictly 
positive-real. There exists e € K, e > and a storage 
function Qq, for the system described in image form 
col(R\,eMR 2 ) with a transfer function (eMR 2 )R~^ 1 , 
and such that Qq, satisfies condition C3 of Theorem\J\for 
the standard SLDS with behaviors 05 1 := ker Ri (^) 
and 05 2 := ker R 2 (^). 

Proof: We begin by computing the state-space 
realization of (MR^R^ 1 corresponding to the stan- 
dard state map (see Lemmas Q] @] and [5J X\ £ 

M niXw [£] for OSi of the form := 

where X 2 £ M( n i- w ) Xw [£] is a state map for 05 2 , and 
X[ € R( ?tl - n2 ) xw [£]. Use the state equation for 05i, i.e. 
compute £Xi(0modiii, to conclude the existence of 
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rii xrii 



and Fx G 



pni xw 



such that AiXi(€) 
-Fi-Ri(f)- Finally, note that since {MR 2 )R\ 



At el 

is strictly proper, we can write M(£)i?2(£) = C\Xi(^) 
for some C\ G M wXni . In the following we call C\ 
the coefficient matrix of MR2 with respect to X\. It 
is readily verified that (A\ , F\ , Ci , 0) is the desired 
realization of (MR 2 )R^ • 

Let e > 0; since (MR^R^ 1 is strictly positive-real, 
also e{MR 2 )R^ 1 is strictly positive-real. Thus, there 
exists aP = P T e Rmx*i such that 



-eCi + F?P 







> 



(13) 



Observe that since the dissipation function for 
col (Pi, eMR 2 ) is a quadratic function of the state X\ 
(see Lemma [3), it holds that PF\ - eCj = 0. Now 
consider the largest solution P + to this LMI, associated 
with (MR 2 )R^ , corresponding to the anti-Hurwitz 
factorization of the polynomial matrix <!?(—£,£) — 
Pi(-£) T M(£)P 2 (£) +R 2 (-0 T M(-0 T Ri(0 where 
= A(— £,,£,). Denote the coefficient matrix of 
the anti-Hurwitz spectral factor with A, i.e. A(^,^) = 
^(C) T ^4(??) with A(£) := AX^). It follows that 
AjP + + P + Ai = -A T A; and from the fact that A\ is 
Hurwitz and A corresponds to the anti-Hurwitz spectral 
factor, that the pair (A\,A) is observable. Now apply 
Theorem 7.4.7 of El to conclude that P + > 0. 

We now show that there exists e such that among the 
solutions P of ( fT"3T ), there exists one with the structure 

Pu -n T p 22 ] 
-P22R P22 



p 



(14) 



where P n G M. n2Xn \ P 22 G R(m-n»)x(»i -"=>), and n G 
g(ni-n 2 )xn 2 represents the gluing conditions involved in 
the switch from ker Pi (4) to ker R 2 (4f), as in Def. 

SI 

Solving (fT3l is a convex problem (see the discussion in 
section 6.3.3 of CQ, the last bullet on p. 93, and ibidem, 
section 2.2.2, line 7 p. 10). Now define the universum 
set U by 



r iXni I P > and P < P + } . 

(15) 



U :={P = P T £ 



Lemma 9. U defined by ([75]) is closed and bounded. 

Proof: See the Appendix. 
Define the feasible set L by 

L:={PeU\3s>0 s.t. PFi - eCj = 0, 

CDfl> holds, P 2 2 > and Pn - n T P 2 2n > j 



(16) 



Lemma 10. The feasible set L defined by rii6D is 
nonempty. 

Proof: See Appendix. ■ 

The results of Lemmas [9] and [10] show that the 
conditions of Assumption 1 p. 43 of Q are satisfied. 
Applying Theorem 16 of Q we conclude that the 
convex optimization problem corresponding to the LMI 
( fT3l ) under the constraints defining L has a solution. 
Now define V((,r]) := ^i(C) 7 PX\(rj); because of 
the structure (fT4l of P, Qqj satisfies condition C3 of 
Theorem [T] ■ 

Since Qq, is a storage function, there exists 
D G IR ,XW [£] such that eR 1 (Q r M(ri)R 2 (ri) + 
sR 2 (C) T M (C) T Pi(r?) = (C + riMC, V) + D(() T D(n). 
It is readily verified that satisfies the conditions CI 
and C2 in Theorem [T] Moreover, in a switch from 53 2 to 
*Bi the value of Qq, remains the same before and after 
the switch, while from Lemma [8] we conclude that in a 
switch from *Bi to Q3 2 it decreases. Qq, is the desired 
Lyapunov function. ■ 

Remark 8. Numerical examples demonstrate that given 
a positive -real completion, the Lyapunov function of 
Theorem |2] does not correspond to a square spectral 
factorization, as was the case for the Lyapunov function 
in Theorem [T] Indeed, the largest storage function for 
MR 2 R^ 1 does not in general have the structure (fT4l 
guaranteeing that its value does not increase at switching 
times. Thus in general a purely algebraic computation of 
the Lyapunov function is unfeasible. However, the proof 
of ((2) (1)) of Theorem |2] shows how the problem 

of computing a Lyapunov function from a given positive- 
real completion can be reduced to solving a structured 
LMI using standard LMI solvers. □ 

Remark 9. There exist pairs of Hurwitz matrices which 
have no positive -real completion, for example the 1 x 
1 matrices: n(£) := 2523677 + 435616£ + 81559£ 2 + 
7000C 3 + 603£ 4 + 24£ 5 + £ 6 and r 2 (f) := 65 + 46£ + 
26£ 2 + 6£ 3 + C 4 - □ 

VI. Conclusions 

We have presented a framework for the analysis of 
linear switched systems in which the different dynamical 
regimes are not described in state-space form, and do not 
share a common state space; pivotal in our approach is 
the concept of gluing conditions, that impose concatena- 
tion constraints on the trajectories of the system at the 
switching instants. We have considered the problem of 
stability in this framework, and we have given Lyapunov 
conditions based on the concept of positive -realness, 
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and constructive methods based on LMIs to compute 
Lyapunov functions. 

The results presented here can only be considered 
preliminary to an application of our approach to other 
problems relevant for switched systems. Some issues 
under study are the establishment of sufficient conditions 
for the stability of switched linear differential systems 
with more than two dynamical regimes; a theory of 
dissipativity; stabilization, and other control problems. 



Appendix 
Proof of technical results 

Proof of Lemma |2]- Let / G £2. Observe that 

/i?2 1 • ^^r 1 = f R x x - Since both f R 2 l and ^^r 1 

are strictly proper, so is their product. Consequently, 

/eli. ■ 

Proof of Lemma \3} To prove the claim on 
the dissipation function Qa> write A(— = 
Q(~0 T Q(0 = D(-0 T N(0 + N(-0 T D(0, and 
define G(£) := A^)^)" 1 and H(£) := QffiDfc)- 1 . 
Then H(-^) T H(^) = G(£) + G(-£) T - Now write 
H(0 =: [M£) T ••• K(C) T ] with h k G M wXl [£], 
k = l,...,w, and observe that because of the strict 
properness of ND -1 it holds that 



lim H(-ju) ' H{ju) 

Ul— >0O 



lim 

LO— >0O 



hx(-ju}) 



[hi(ju) ■■■ h r (jLo)} 



h r (-ju) T _ 
= lim G{-juj) T + GCju) = . 

The diagonal entries of H(—ju)) T H(juj) 
are ||/ii(jw) || 2 , ... , \\h v {juj)\\ 2 , and thus 
lim aJ _ s . 0O H/ifc^'a;)!! 2 = for k = 1, ...,w. It follows that 
the columns h k (Q of H(£) = QfyDfe)* 1 are strictly 
proper, as was to be proved. 

To prove the statement about the storage func- 
tion Qqj with associated dissipation rate Q/a> write 

Q(0 T Q(ri) + (C + vMCv) = d(c) t n(v) + 

N(() T D(r]). Premultiply both sides by D(()~ T to ob- 
tain D(C)- T Q(0 T Q(r]) + (C + v)D(0~ r *(C,V) = 
N(r]) + D(()- T N(() T D(r]). Now compute 

lim D(Cy T Q(0 T Q(v) + (C + v)D(cr T *(c, v) 

= lim N{rj) + D(()~ T N(() T D(r]) . 

Since ND^ 1 and QD -1 are strictly proper we obtain 

lim c ^ 00 (C + ??) J D(C)~ T ^(C,??) = lim^ooiV^). Using 



the given canonical factorization of ^(£,77) we obtain 
lim J D(C)" T A/(C) T (CS* + r ? S^)G(r ? ) 

s-oo 

= lim TV (77) = TV (77) . 

s-oo 

Assume by contradiction that i?(C) := ^(C)-C(C) -1 is 
not strictly proper, and write H(Q = Hl( l + ... + ifo + 
S(C), where L > and 5(C) is strictly proper. Then 
lim^oo D(()- T G(() T (££* + 7/S*)G(r7) is infinite, a 
contradiction. Consequently, := M(£)-D(£) _1 is 

strictly proper. ■ 
Proof of Lemma |U - The claim on the limit being 
finite follows from the fact that X^R^ 1 is strictly proper. 
Recall from section IlLCTl that there exist A\ G JR niXn \ 
F 1 G M™i>< w such that 



(17) 



Now observe that since R2R1 1 is strictly proper and i?i 
is column-reduced, each of the column degrees of R2 
is at least one less than the degree of the corresponding 
column of R\. It follows that lim^oo £X2(£)Ri(£) = 
O n2 xw ; and consequently the first n<i rows of F\ are equal 
to zero, i.e. 



lim 








"0" 


cximiicr 1 . 




k 



(18) 



with k G rK-™2)xw_ 

Since R2RT 1 is strictly positive real and using the 
equalities *(C,r?) = Xi(C) T $Xi(r/) and A(C,r?) = 
Xi(C) T AXi(r/) established previously, write 

(C + 7 7 )X 1 (C) T *X 1 (r 7 ) + Xi(C) T AX 1 (r ? ) 
= (AiXifO + Fii?!^)) 1 "*^!^) 

+ Xi(0 T * (A 1 X 1 { V ) + FxR^r,)) + X 1 (C) t AX 1 (t 7 ) 
= X!(C) T + X x (r/) + i?i(C) T J Pi T *X 1 (r ? ) 

+ Xi(C) T *Fii2i(7/) + Xi(C) T AJTi(7/) 
= i? 1 (C) T i?2W +i? 2 (C) T i?i(r/) • 

Now multiply on the left by i?i(£)~ T both sides of the 
equality, and take the limit for £ — > 00; since XiR^ 1 
and R2R1 1 are strictly proper, we obtain i?2(0 = 

FJ^x^e) = & T (*7 2 x 2 (0 + *22^(0)- ■ 

Proof of Lemma^s Following the same argument as 
in the proof of Lemma |4] we write A\X\{£) = £Xi(£) — 

Ri(£). From Lemma|4]and the fact that R2 is square 

and nonsingular, it follows that k T is of full row rank, 
and consequently n\ — 712 > w. 

In order to prove that n\ — 122 = w, assume by 
contradiction that n\ — 112 > w. Then there exists one 
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column of Ri whose column degree is greater than or 
equal to the degree of the corresponding column of R2 
plus two; w.l.o.g. assume that this is the first column of 
R\. Then, possibly rearranging the rows of X[, write 



St-l 



Moreover, permute the rows of X[ so that e^* 5 * -1 , 
i = 1, . . . , w, are its last w rows. Note that these 
operations correspond to permutations of the rows of 

fc = iim^ 00 ex((e) J Ri(6~ 1 - 

Given this structure of X[, note that the last w 
rows of k contain the inverse of the highest column 
coefficient matrix of R±, while its first n\ — 712 — w 



rows are equal to zero, i.e. k 



("I 



-n 2 — w)Xw 



k' 



with k' G ]R wXw nonsingular. Denote the last w rows 
of ^ 2 2 by [% 2 ¥£ 2 ] with % 2 G R wX K-«2- w ) 
and G JR wXw . From Lemma @J P 2 (£) = 

k T (*J 2 X 2 (t) + # 22 X((£)), and consequently R 2 (£) = 
k' T ^ 2 X 2 (0 + £;' T [% 2 tfg 2 ] X[(0, where e 
R wxn2 denotes the last w rows of ^>J 2 . Note tnat due 
to the positive defmiteness of $22, also ^ 22 > 0; then 
the highest column coefficient matrix of P 2 (£) is k'^ 22 , 
and its first column is associated with the column degree 
8\ — 1 = (5^ + 1. This contradicts the fact that the 
first column of P 2 (£) has degree 5\. Consequently we 
conclude that n\ — 17,2 = w and that k is nonsingular. The 
claim is proved. ■ 
Proof of Lemma® From the definition of IT it fol- 



lows that [IT 



-n 2 )] 



X((£)modP 2 



X((0. Write X((0i? 2 (0^ = ^(0 + S(0, with 
S'(^) a strictly proper polynomial matrix and P G 

RCni-naJxw^]. then 



We now compute ((^ + rj)^ ((, rj) mod Ri =: — A(£, 77). 
Since $22 is by assumption positive definite and sym- 
metric, there exists a nonsingular matrix G such that 
$ 22 = G T G. Thus 

(C + 7?)*(C,7?) =(C + r?)* / (C,7?) 

+ Ci? 2 (C) T ^(C) T G' T GP( ?? )^2( ? ?) 
+ P 2 (C) T P(C) T G T GP(r ? )P 2 (r ? )7 ? . 

$'(C,^) is a function of the state of *B 2 , i.e. ^'(C) 7 ?) = 
X 2 {C) T ^' X 2 {r]) for some $' G M n2Xri2 , and ^^r 1 
is strictly proper. Consequently, in order to com- 
pute (C + r])^((, rj) mod Pi it is enough to com- 
pute GP(0^2(0? mod ^i- From Gil) it follows that 
GP(0R2(m = G[-U I {ni _ n2) ] X^; moreover 
by 03, = MX^) + Fifli(e). Thus 

gp(£)p 2 (6£ =g [-n i (ni _ re2) ] (^^(e) + FiiMC)) 
=G[-n / (ni _„ 2) ] AiXiCO 
+ G[-n j^^jjFiUiCO 



[n 



( (ni-n 2 ) 



jcj(0-P(0fl2(0 



(19) 



Substituting OH) in QUI) , we obtain OB- ■ 
Proof of Lemma 0" From (fTTI ) it follows that 

(C + T7)*(C, 77) =(C + r?)^(C,r?) 

+ (C + 7 ? )P2(C) T ^(C) T ^22P(^)P 2 (^) 

=(C + ^)*'(C^) 

+ CP2(C) T ^(C) T ^22P(r?) J R2(r?) 

+ P 2 (C) T J P(C) T *22P(??)P2(??>?. 



--G[-U J, 



(tii- n 2 ) 



A\X\(£) + GBRi(£) . 



where 5 G rK-" 2 ) Xw Then 
GPiOMOC^odRx = G [-11 / (ni _ n2) ] AiXi(0 

Note that P satisfies B = lin^-wo P(£)P 2 (£)£Pi(£) _1 . 
Consequently, we obtain 

-A(C,t?) 

= (C + r ? )$'(C,r ? ) 

+ (CP 2 (C) T P(C) T - J Ri(C) T 5 T ) G T GP(r,)R 2 (r,) 

+ R 2 (C) T P(C) T G T G (P( V )R 2 ( V )r) - BRM) 

= (C + r/)^(C, r/) + CR2(C) T P(C) T G T GP(r ] )R 2 (7 ] ) 

-R 1 (C) T B T G T GP( V )R 2 (v) 

+ R2(0 T P(0 T G t GP(7 1 )R2(t 1 )v 
-R 2 (C) T P(0 T G T GBR 1 ( v ) 

= (C + »7)*'(C, »7) + (C + v)R2(0 T P(C) T G T GP(r ] )R 2 (r ] ) 
-R l {C) T B T G T GP{i 1 )R 2 {i 1 ) 
-R2(() T P(() T G T GBR 1 (7 ] ) . 

From OB we conclude that (C + r/)$ / (C,r/) + (C + 
V )R 2 (C) T P(0 T G T GP( V )R2(v) = (C + vMCvl and 
consequently 

-A(C,77) =(C + r])*(C + V) 

-Pi(C) T P T ^2 2 P(??)P2(r?) 

-P 2 (C) T J P(C) T *22PPl(??) , 

as was to be proved. ■ 
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Proof of Lemma® To prove closedness, fix a basis 

{Pi € R niXni , Pi = Pj} i=1 ni(n 1+ i) for the space of 

2 -p 
symmetric ni -dimensional matrices, i.e. every P = P 

can be written as P = Yl<i=x 2 a iPi f° r some coeffi- 
cients aj € R. P G J7 if and only if 



"l("l+D 



Oni xni 
Ojii xrti 



Om xni 

p+ 



+ 



CY; 



P 







0, 



rti xni 
Pi 



The left hand side of 



'ni xni 

(20) 

is a matrix whose entries 



are real linear polynomials in the indetermi- 
nates aj. The set consisting of all the rai ^" 1+1 ^ -ples 
cki, . . . , Q ni(ni+i) such that d20l holds is a spectrahedron, 
and it is a closed set, see p. 1944 of JSJ. Consequently 
?7 is closed. 

To prove boundedness, choose the metric of the in- 
duced Euclidean norm on the space of matrices; then U 
is contained in the closed ball with radius equal to the 
largest eigenvalue of P + > 0. 

Proof of Lemma [70|' Partition C\ =: Nx N 2 

with Nx £ R wXn2 and iV 2 € R^K-"*), 



i.e. 



CxXx(0 = 



Define k e r("i by pT = . fe Tj. then 

P of the form dHJ satisfies PP X - Cj = if and 

only if -n T P 22 A: -^7 = ° and ^22 fc - JVj" = 0. 
We now show that if the second equality holds, then 
the first equality is automatically satisfied. Rewrite equa- 
tion (ED as N 2 X[{£) = M(0%(0 - NxX 2 (0 and 
compute iV 2 X((£)modP 2 = -NxX 2 (£). Recall that 
X((£)modP 2 = ILY 2 (£), and conclude that 



iVi N 2 



MO 



(21) 



n T iv 2 T 



N 



(22) 

Consequently if P 22 A; = Nj , then (1221 is automatically 
satisfied. 

Thus to prove that L is non-empty it suffices to show 
that there exists e > such that among the matrices P 
with structure (fT4l . and such that k T P 22 = eN 2 , there is 
at least one such that also P 22 > 0, Pq - n T P 22 n > 0, 
and that is smaller than or equal to the maximal storage 
function P + . 

In order to do this, transform congruently P in (fT4l 

. , \Pxx - n T p 22 n o " 

to the form 

U p2 2 _ 

transformation to P + . To keep the notation light, in 
the following we use the same symbols to denote the 
matrices P and P + before and after the change of 
variables. Partition the (congruence-transformed) P + as 



, and apply the same 



P+ 



where P^ S 



xn 2 



p+ p+l 
r X\ r X2 

p+T p + 
r X2 r 22. 

{n 1 -n 2 )x(n 1 -n 2 ) p+ ^ ^n 2 x(n 1 -n 2 ) 

Let e be a positive real parameter and define 



22 



p22,£ 
Pll, £ 



eP 



22 



(1 + e)n T P 22je n = (1 + e)en T P+n(23) 



> q Note that since P 2 + 2 > 0, also Pxx, E > and P 22j£ > 0. 
We now prove that there exists £ £ R, e > 0, such that 



r p + 


P+~ 

r X2 




"Pll, £ 


p+T 
+ 12 


P+ 
^22. 







n T p 22 , £ n 



p + 
r xx 



e 2 n T p+n 

p+T 
r X2 



(1 



12 





p22,£ 
> . 



so, compute factorizations P 22 = V T V, 
W T W for some nonsingular W, V of 
dimension nx — n 2 and n 2 , respectively. Rewrite the last 
inequality as 



To do 

n T p+n 



p+T 
V12 



12 
^+ 
22 





" 












e 2 /, 



"2 







eP 



TT" 





> 



then bring it to the equivalent form 

\W~ T 






v-' 



e 2 l 




n-2 





r p+ p+- 

Ml ^12 




W~ x 




p+T p + 
+ 12 r 22. 




o y- 1 







> o . 




ll- 


-n 2 _ 







The first matrix on the LHS is known, and it is positive- 
definite; call it P +/ . To satisfy the inequality, choose < 
£ < A min (P +/ ), the latter being the smallest eigenvalue 
of P +/ . 

In order to complete the proof of the Lemma, we 
need to show that having chosen e in this way and 
defined P from (l23l l. we obtain a matrix satisfying the 
equality PF\ — eCj = 0. Observe that P + before the 
congruence transformation satisfied P + Pi — Cj = 0, 
and consequently before the congruence transformation 
k 1 P 22 = Nj; it follows from the definition of P 22j£ in 
(123T ) that k T P 22)£ = eN 2 . Using <[22]> we conclude that 
PPi - eCj 



0. 



References 

[1] S. Boyd, L.E. Ghaoui, and V. Balakrishnan. Linear Matrix In- 
equalities in System and Control Theory. Society for Industrial 
Mathematics, 1994. 

[2] W.A. Coppel. Linear Systems, Notes in Pure Mathematics, 
Australian National University. 6, 1972. 

[3] J. P. Hespanha and A. S. Morse. Switching between stabilizing 
controllers. Automatica, 38(11): 1905-1917, 2002. 

[4] T. Kailath. Linear Systems. Prentice-Hall Information and 
System Sciences Series, 1980. 



14 



[5] D. Liberzon. Switching in Systems and Control. Birkhauser. 
Boston, Basel, Berlin, 2003. 

[6] H. Lin and P.J. Antsaklis. Stability and stabilizability of 
switched linear systems: A survey of recent results. IEEE 
Transactions on Automatic Control, 54(2):308-322, 2009. 

[7] A. Nedic. Lecture notes on Optimization I. Network Mathe- 
matics Graduate Programme Hamilton Institute. 

[8] T. Netzer, D. Plaumann, and M. Schweighofe. Exposed faces 
of semidefinitely representable sets. SIAM J. Control Optim., 
20(4): 1944-1955, 2010. 

[9] R. Newcomb. Linear Multiport Synthesis. New York: McGraw- 
Hill, 1966. 

[10] R. Peeters and P. Rapisarda. A two-variable approach to solve 
the polynomial Lyapunov equation. Syst. Contr. Lett., 42:117- 
126, 2001. 

[1 1] J.W. Polderman and J.C. Willems. Introduction to Mathematical 
System Theory: A Behavioral Approach. Springer, Berlin, 1997. 

[12] P. Rapisarda and P. Rocha. Positive realness and Lyapunov 
functions. Proceedings of the 20th International Symposium 
on Mathematical Theory of Networks and Systems, 2012. Mel- 
bourne, Australia. 

[13] P. Rapisarda and J.C. Willems. State maps for linear systems. 
SIAM J. Control Optim., 35(3): 1053-1091, 1997. 

[14] P. Rocha, J.C. Willems, P. Rapisarda, and D. Napp. On the 
stability of switched behavioral systems. 50th IEEE Conference 
on Decision and Control and European Control Conference 
(CDC-ECC), pages 1534-1538, 2011. 

[15] R. Shorten, K.S. Narendra, and O. Mason. A result on common 
quadratic Lyapunov functions. IEEE Transactions on Automatic 
Control, 48(1): 110-1 13, 2003. 

[16] R.N. Shorten and K.S. Narendra. On common quadratic 
Lyapunov functions for pairs of stable lti systems whose sys- 
tem matrices are in companion form. IEEE Transactions on 
Automatic Control, 48(4): 6 18-621, 2003. 

[17] Z. Sun and S.S. Ge. Switched Linear Systems: Control and 
Design. Springer- Verlag, New York, 2005. 

[18] S. Trenn. Switched differential algebraic equations. Dynamics 
and Control of Switched Electronic Systems. Chapter 6 of: 
Francesco Vasca and Luigi Iannelli (eds.), Springer Verlag, 
2012. 

[19] H. L. Trentelman and P. Rapisarda. New algorithms for 
polynomial j-spectral factorization. Math. Contr. Sign. Syst., 
12:24-614, 1999. 

[20] J.C. Willems and H.L. Trentelman. On quadratic differential 
forms. SIAM J. Control Optim., 36:1703-1749, 1998. 

[21] H. Ye, A.N. Michel, and L. Hou. Stability theory for hybrid 
dynamical systems. IEEE Transactions on Automatic Control, 
43(4):461^174, 1998. 



